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We discuss some aspects of a recently proposed semi-classical transport theory 
for QCD plasmas based on coloured point particles. This includes the derivation 
of effective transport equations for mean fields and fluctuations which relies on 
the Gibbs ensemble average. Correlators of fluctuations are interpreted as colli- 
sion integrals for the effective Boltzmann equation. The approach yields a recipe 
to integrate-out fluctuations. Systematic approximations (first moment, second 
moment, polarisation approximation) based on a small plasma parameter are dis- 
cussed as well. Finally, the application to a hot non-Abelian plasma close to ther- 
mal equilibrium is considered and the consistency with the fluctuation-dissipation 
theorem established. 

1 Introduction 

Reliable theoretical tools to describe transport phenomena of hot or dense 
QCD plasmas in- or out-of equilibrium are at the basis for an understanding 
of the quark-gluon-plasma. This extreme state of matter where quarks and 
gluons are no longer confined is expected to be produced within the up-coming 
heavy ion experiments at RHIC and LHC. 

In the present contribution a recently proposed semi-classical method to 
derive effective transport equations for QCD is reviewed, developed in collab- 
oration with Cristina Manuel in_ref. □□. The interest in a formalism based 
on a classical transport theory H is that main properties of a hot quantum 
plasma can already be understood within simple classical terms. Indeed, the 
soft non-Abelian gauge fields - those having a huge occupation number - can 
be treated as classical fields, while the hard gauge field modes and the quarks 
can be treated as classical particles. This approach has been known since long 
for Abelian plasmas a, but a consistent extension to the non-Abelian case was 
longtime missing. 

Here, we discuss some conceptual aspects of this method, and in particular 
the step from a microscopic to a macroscopic (kinetic) formulation, the close 
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relationship between the ensemble average and the notion of coarse-graining, 
and systematic expansions in a small plasma parameter or gauge coupling. The 
latter is at the basis for the integrating-out of fluctuations and the derivation of 
collision integrals. For a hot plasma close to equilibrium it has been shown that 
the HTL effective theory follows to leading^ order in an expansion in moments!!!. 
We review how Bodeker's effective theoryB is recovered at second order in this 
expansion. Finally, the compatibility with the fluctuation-dissipation theorem 
is established^. 

2 Classical particles and fields 

The starting point for a classical transport theory is to consider an ensemble of 
classical point particles interacting through classical non-Abelian gauge fields. 
Such an approach is known since longJbr Coulomb plasmas of charged point 
particles interacting through photons □. The new ingredience in the case of 
QCD is that the particles have to carry a non-Abelian colour charge Q a , where 
the colour index runs from a = 1 to TV 2 — 1 for a SU(7V) gauge group. These 
particles interact self-consistently amongst each others, that is, through the 
classical gauge fields created by the ^articles. Their classical equations of 
motion have first been given by Wongla, 

m=- = P* , m-f- = gQ a F^p v , m-f- = -gf abc p»A»O c . (1) 
or dr dr M 

Here, denotes the gauge field, F* v = d^A" - d u A^ + gf ahc A\A C V the corre- 
sponding field strength and z = (x,p,Q)(t) the world line of a particle. The 
particles interact self-consistently through classical Yang-Mills fields, 

= ju ^ (2) 

where the current J is the sum of the currents for the individual particles, 
J = gJ2i I d,TS(x — Xi)dxi/dr. 

To further a kinetic description it is useful to introduce a microscopic one- 
particle distribution function f(z) ~ J drS[z — Zi(r)] which describes an 
ensemble of such particles. Making use of (Q), its kinetic equation iscl 

V» - 9f*°AlQrJL _ gQaKv f(x , p , Q) = . (3) 

This Boltzmann equation is collisionless, since df /dr = (Liouville's theorem). 
However, it contains effectively collisions due to the long range interactions as 
shall become clear in the sequel. The Wong equations transform covariantly 
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under gauge transformations which entails that / and (Q) are gauge invariant!]. 
Eq. (||) is completed by the Yang-Mills equation (Q), where the current J 
becomes now a functional J[f] of the one-particle distribution function. These 
equations constitute the basis for the construction of a semi-classical kinetic 
theory. 

3 Kinetic theory 

Within the semi-classical approach all informations about properties of the 
QCD plasma are given by the microscopic dynamical equations (||) and (||). 
However, for most situations not all microscopic informations are of relevance. 
Of main physical interest are the characteristics of the QCD plasma at macro- 
scopic scales. This includes quantities like damping rates, colour conductivities 
or screening lengths within the kinetic regime, or transport coefficients like 
shear or bulk viscosities within the hydrodynamic regime. The microscopic 
length scales, like typical inter-particle distances, are much smaller than such 
macroscopic scales. 

Given the statistical ensemble - representing the state of the system - 
its macroscopic properties follow as functions of the fundamental parameters 
and the interactions between the particles. Within kinetic theory, the basic 
'macroscopic' quantity is the one-particle distribution function from which all 
further macroscopic observables can be derived. The aim of a kinetic theory is 
to construct - with as little restrictions or assumptions as possible - a closed 
set of transport equations for this distribution function. Such an approach 
assumes implicitly that the 'medium', described by the distribution function, 
is continuous. If the medium is not continuous, stochastic fluctuations due 
to the particles should be taken into account as well, and their consistent 
inclusion leads to a coupled set of effective transport equations for correlators 
of fluctuations and the one-particle distribution function u. In this classical 
picture, / is a deterministic quantity once all initial conditions have been fixed. 
Unfortunately, this is not possible, and / has to be considered as a stochastic 
(fluctuating) quantity instead. The random fluctuations of the distribution 
function are at the root of the dissipative character of the effective transport 
theory. 

4 Coarse graining and ensemble average 

The step from a microscopic to a macroscopic (kinetic) description requires an 
appropriate definition of macroscopic quantities in the first place. There are 
(at least) two possible options. The first one consists in taking both volume 
and time averages of the microscopic distribution function /, the non-Abelian 
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fields and thp Boltzmann equation (j3|) over characteristic physical length and 
time scales 0. Clearly, the notion of a characteristic physical scale depends 
on the particular problem studied. A reasonable choice for a coarse graining 
scale r p h has to met some few criteria. The coarse graining scale r p h should 
be sufficiently large such that the number of particles -/V p h contained in r^ h is 
3> 1. This is at the basis for a plasma description in the first place, because 
N p h 3> 1 ensures that many particles interact coherently within a coarse- 
graining volume (quasi-particle behaviour). The inverse of iV p h is also known 
as the plasma parameter e. In addition, r p h should be larger than a typical 
two-particle correlation length r corr . This ensures that two-particle correlators 
<72 can be neglected. Finally, r p h should be smaller than the macroscopic scales 
of the plasma to be investigated, like typical relaxation lengths r ro i. This way 
it is guaranteed that the scales of interest are not washed-out. Such coarse 
grainings remove (irrelevant) microscopic information and thereby modify the 
transport equation which is expected to become of the Boltzmann-Langevin 
type, 

P» (D„ - gQaF^d;) f = C[f] + C ■ (4) 

Here, we introduced Dfj\A\f = [<9 M — gf^QcA^.bd®}/ and the shorthands 
dfj, = d/dx 11 , = d/dp^ and = d/dQ a . The new terms on the r.h.s. are 
the collision integral C[f] and a related source for stochastic noise £ on the 
r.h.s. of (^) are due to the fact that the Boltzmann equation (|^) is quadratic 
and cubic in the fields. Physically speaking, these terms arise because the 
coarse-graining integrates-out the short range modes within a coarse-graining 
volume. The interactions of such modes can result into additional effective 
interactions for the remaining long range modes. 

A second route for obtaining kinetic equations consists in taking an en- 
semble average of the microscopic transport equation. In the present con- 
text we consider particles in a phase space, hence the appropriate average is 
the Gibbs ensemble average henceforth denoted by (•••). All statistical infor- 
mations of the system are then contained in the correlators {5f ■ ■ ■ Sf) with 
8f = f — (/} = /— / denoting the statistical fluctuations of / about its mean 
value. The most important equal-time correlator is the quadratic one given by 

(*/x 1 p,o«5/xV,<y) = (2vr)V 3 )(x - x')<5 (3) (p - p')S(Q - Q') f + 0(g 2 ) (5) 

where we have neglected higher order corrections due to the two-particle corre- 
lation function gi . Similar formula hold for higher order equal-time correlation 
function. Notice that (|^) does not require further informations about the state 
of the system and applies equally to systems in- or out-of equilibrium. While 
has been proven rigorously for particles with classical statistics, its gener- 
alisation to the case of quantum statistics is known only for weakly interacting 
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particles and consists in replacing 

/-»/(l±/) (6) 
in (JsJ) for Bose-Einstein (+) or Fermi-Dirac (— ) statistics. 

5 Effective transport equations 

We now apply the Gibbs ensemble average to the Boltzmann equation (||) in 
order to derive an effective kinetic equation of the form (|]). To that end, we 
split / into its mean value and a fluctuating part 

f(x,p,Q) = f(x,p, Q) + 5f(x, p, Q) , (7) 

where / = (/}. The mean value of the statistical fluctuations vanish by def- 
inition, (Sf) = 0. This separation into mean fields and statistical, random 
fluctuations corresponds effectively to a split into long vs. short wavelength 
modes associated to the mean fields vs. fluctuations. In addition to (0) we 
split the gauge field into a mean field part and fluctuations, 

A»(x) = A^{x)+a' i (x) , (8) 

with A = (A) and (a) = 0. This split should been seen on a different footing 
as the split (0) because the Gibbs ensemble average is defined for the particle- 
like degrees of freedom in the first place. One may wonder whether (|J) - and 
the subsequent split of the dynamical equations - is compatible with gauge 
symmetry because the gauge field transforms non-lincarly. The situation is 
very similar to the background field method used in path integral formulations 
of QFT and the compatibility of (||) n with gauge symmetry follows along the 
same lines. For more details, see ref. 

Coming back to the initial set of dynamical equations (|^) and (|^), we 
perform the ensemble average to find 

V» - gQaW) / = (r?> + (0 ■ (9) 

for the effective transport equation, and 

D ll F t " / + (J»J = r . (10) 

for the effective Yang-Mills equations. In (^|) and ([l0]), we collected all terms 
quadratic or cubic in the fluctuations into the functions 

n = gQaV' 1 {iP^ v - D v a„) a + gr^al] d%Sf(x,p,Q) , (11) 
£ = gp»f abc Q c [d^al Sf(x,p, Q) + ga«ald;f(x : p, Q)] , (12) 
JZ = g [f dbc D^ ad a b ^ c + / a6c a 6;M - W l ) c + <7/ cde «)] • (13) 
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The set of equations (^) - (|T3|) is not yet closed as they still involve correlators 
of fluctuations. Hence, we need their dynamical equations as well. Subtracting 
(||) from (||) we find for the dynamics of Sf 

V 11 {p„ - gQ a F« v d») Sf = gQ a (D„a v - D^p^J 

+9P"a b4i r bc Q c dQf + r) + £ - (r, + . (14) 

Subtracting (||) from (|l^) yields the corresponding equation for the gauge field 
fluctuations, to wit 

(5V - D»{D„a v )) a + 2gf abc Fg v a c , v + J££ - (J££) = SJ a ^ . (15) 

Eqs. ( |l4] ) and (|l^) are the master equations for all higher order correlation 
functions of fluctuations, e.g. multiplying ( jl4| ) by Sf and taking the ensem- 
ble average yields the Boltzmann equation for the correlator (SfSf). This 
hierarchy of dynamical equations is similar to the BBGKY hierarchy. All the 
equal-time correlators of Sf can be derived from the basic definition, of the 
Gibbs ensemble average in phase space and serve as initial conditionsH. 

We remark that the set of dynamical equations @ - (|l5|) is equivalent 
to the original set of microscopic equations. It is exact in that no further 
approximations - apart from the assumptions leading to ([!]) and (||) - have 
been performed so far. Furthermore, it can be applied to out-of-equilibrium 
situations simply because the Gibbs ensemble average is not bound to a plasma 
close-to-equilibrium. 

6 Collision integrals 

Let us comment on the terms r], £ and J nnc in the dynamical equations. In 
the effective Boltzmann equation, the functions (rj) and (£) appear only after 
the splitting (Q) and (||) has been performed. These terms are qualitatively 
different from those already present in the transport equation. The correlators 
{rf) and (£) are interpreted as effective collision integrals of the macroscopic 
Boltzmann equation. The fluctuations in the distribution function of the par- 
ticles induce fluctuations in the gauge fields, while the gauge field fluctuations, 
in turn, induce fluctuations in the motion of the quasi-particles. In the present 
formalism, the correlators of statistical fluctuations have the same effect as 
collisions. This yields a precise recipe for obtaining collision integrals within 
semi-classical transport theory. 

In the Abelian limit the transport equations reduce to the known set of ki- 
netic equations for Abelian plasmas and only the collision integral (77) survives. 
Here it is known that (rj) can be explicitly expressed as the Balescu-Lenard 
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collision integral!]. This substantiates the correspondence between fluctuations 
and collisions in an Abelian plasma. 

At the same time we observe the presence of stochastic noise in the effective 
equations. The noise originates in the source fluctuations of the particle dis- 
tributions and induces field-independent fluctuations to the gauge fields. The 
corresponding terms in the effective Boltzmann equations are therefore rj, £ 
and J fluc at vanishing mean field or mean current. 

Finally, we observe the presence of a fluctuation-induced current ( J fluc ) in 
the effective Yang-Mills equation for the mean fields. This current, due to its 
very nature, stems from the induced correlations of gauge field fluctuations, ft 
vanishes identically in the Abelian case. While the collision integrals are linear 
in the quasi-particle fluctuations, the induced current only contains the gauge 
field fluctuations. As the fluctuations of the one-particle distribution function 
are the basic source for fluctuations, we expect that a non-vanishing induced 
current will appear as a subleading effect. 

In order to find explicitly the collision integrals, noise sources or the 
fluctuation-induced currents for non-Abelian plasmas, one has to solve first 
the dynamical equations for the fluctuations in the background of the mean 
fields. This step amounts to 'integrating-out' fluctuations. In general, this is 
a difficult task, in particular due to the non-linear terms present in 77, £ and 
J fluc . This will only be possible within some approximations. 

7 Plasma parameter and moment expansions 

When it comes to solving the kinetic equations - and in particular the fluctua- 
tion dynamics - it is necessary to identify a small expansion parameter which 
allows for systematic approximations of the hierarchy of dynamical equations 
for correlation functions. The non-Abelian plasma is characterised by two di- 
mensionless parameters, the gauge coupling g and the plasma parameter e. A 
small plasma parameter is mandatory for a kinetic description to be viable. For 
a classical plasma, e is an independent parameter related to the mean particle 
number and the gauge coupling. For a quantum plasma, the mean particle 
number cannot be fixed arbitrarily, and e scales proportional to the gauge cou- 
pling as e ~ j 3 C 1, explaining why a kinetic description for quantum plasmas 
is intimately linked to the weak coupling limit. 

Here, two systematic approximation schemes are outlined: an expansion 
in moments of the fluctuations and an expansion in a small gauge coupling. 
Although they have distinct origins in the first place they are intimately linked 
due to the requirements of gauge invariance. 

A systematic perturbative expansion in powers of g can be done because 
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the differential operator appearing in the effective Boltzmann equation (£)]) 
admits such an expansion. The force term gp^Q a Fu U ^p i s suppressed by a 
power of g as compared to the leading order term p^D^. In this spirit, we 
expand 

/ = / (0) +5/ (1) + ... (16) 

This is at the .basis for a systematic organisation of the dynamical equations 
in powers of gu To leading order, this concerns in particular the cubic corre- 
lators appearing in (77) and (J fluo ). They are suppressed by a power of g as 
compared to the quadratic ones. At the same time, the quadratic correlator 
~ fabc(o b a c ) within (£) is also suppressed by an additional power of g and 
should be suppressed to leading order. 

An expansion in moments of the fluctuations has its origin in the frame- 
work of kinetic equations, which describe the coherent behaviour of the parti- 
cles within some physically relevant volume, described by the plasma parame- 

—1/2 

ter e. The fluctuations in the number of particles ~ iV ph become arbitrarily 
small if the physical volume - or the number of particles contained in it 
can be made arbitrarily large. For realistic situations, both of them are finite. 
Still, the fluctuations remain at least parametrically small and suppressed by 
the plasma parameter. Hence, the underlying expansion parameter for an 
expansion in moments of fluctuations is a small plasma parameter 

e«l. (17) 

The leading order approximation in an expansion in moments assumes e = 0. 
This is the first moment approximation which consists in imposing 

/ = /, or 6f = 0, (18) 

and neglects fluctuations throughout. Sometimes it is referred to as the mean 
field or Vlasov approximation. It leads to a closed system of equations for the 
mean one-particle distribution function and the gauge fields. In particular, the 
corresponding Boltzmann equation is dissipationless. 

Beyond leading order, the second moment approximation takes into ac- 
count the corrections due to correlators up to quadratic order in the fluctua- 
tions (Sfdf). All higher order correlators like (SfiSfz ■ ■ ■ Sfn) = for n > 2 arc 
neglected within the dynamical equations for the mean fields and the quadratic 
correlators. This approximation is viable if the fluctuations remain sufficiently 
small. We remark that the second moment approximation no longer yields a 
closed system of dynamical equations for the one-particle distribution function 
and quadratic correlators because the initial conditions for the evolution of 
correlators, which are given by the equal-time correlation functions as derived 
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from the Gibbs ensemble average, involve the two-particle correlation functions 
g2- Hence, we have to require that two-particle correlators remain sufficiently 
small as compared to products of one-particle distribution functions. This 
approximation is known as the approximation of second correlation functions, 
sometimes also referred to as the polarisation approximation^. 

For the dynamical equations of the fluctuations these approximations im- 
ply that the terms non-linear in the fluctuations should be neglected to leading 
order, setting 

V-(v)=0, £-<0 = 0, J f .uo - (Jnuo) = . (19) 

It permits truncating the infinite hierarchy of equations for the mean fields and 
the correlators of fluctuations down to a closed system of differential equations 
for both mean quantities and quadratic correlators. The polarisation approx- 
imation is the minimal choice necessary to genuinely describe dissipative pro- 
cesses, because it takes into account the feed-back of stochastic fluctuations 
within the particle distribution function. 



Finally we notice that small fluctuations Sf about / in (16) are already 
of the same order of magnitude as the small corrections from g /W once 
gf^/f^ ~ ^ 1- This means that fluctuations do contribute to the 

effective kinetic equations even if only small departures from the mean value 
.g/*- 1 -* <C are considered. This shall become more explicit at the example 
of a plasma close to thermal equilibrium. 



8 Thermal plasmas 

We now turn to the specific example of a hot plasma close to thermal equi- 
librium. 'Hot' implies that particle masses can be neglected to leading order 
to <C T, and that the gauge coupling, as a function of temperature, is very 
small g(T) < 1. To simplify the analysis we shall perform some approxima- 
tions, all of which can be understood as a systematic expansion in g (and 
e). After ensemble averaging, the distribution function / is effectively coarse- 
grained over a Debye volume. Fluctuations of / within a Debye volume are 
parametrically suppressed by powers of g. We expand / about the equilibrium 
distribution function / oq to leading order in g as 

f = f cq + gf {1) +Sf ■ (20) 

Solving (|]) in the firsLsnoment approximation, that is using ( Bfl ) for Sf = 0, 
has been shown in ref . to reproduce the HTL effective theory. 

Beyond the HTL level we employ the polarisation approximation where 
Sf 7^ 0. It consists in discarding cubic correlator terms in (^) in favour of 
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quadratic ones. Employing ([19|) means that the effects of collisions are ne- 
glected for the dynamics of the fluctuations. All these approximations can be 
systematically improved to higher order. 

Let us consider colour excitations, described by the colour current density 

J£(z,v) = ^ J dQdp d\p\\p\ 3 e( PQ )S(p 2 )Q a f(x,p,Q) , (21) 

and = p^/po — (l,v),v 2 = 1. The current J is obtained after an angle 
average over the directions of v as J[x) = J ^-J(x, v). The colour measure 
dQ is normalised J dQ = 1 and contains the Casimir constraints of the gauge 
group, such as J dQQ a Qb = C^ab with the quadratic group Casimir Ci = N 
for particles in thft-adjoint representation of SU(N) and C*2 = \ for particles in 
the fundamentalH'Q. The approximate dynamical equations for the fluctuations 
become 

8J") a = -m 2 D v»v» (D^a - D aJ a - g f abc v» a\j^ , (22) 
(£V - D»(Da)) a + 2gf abc F£»a c , v = 6J% , (23) 
where the quantum Debye mass for SU (N) with Np quarks and anti-quarks 

2 /-oo / Foq Jeq 

m 2 D = -^ dpp 2 (Nlf + N F lf 
7T Jo \ dp dp 

= g 2 T 2 (2N + N F )/6 (24) 

obtains from the equilibrium (Bose-Einstein or Fermi-Dirac) distribution func- 
tion and the group representation of the particles. The classical Debye mass 
follows by using the Maxwell-Boltzmann distribution, instead. Solving for the 
fluctuations in the present approximations yields a closed expression for S^L. 
and an iterative expansion in powers of the background fields for a (see ref. □ 
for the explicit expressions). The seeked-for dynamical equation for the mean 
quasi-particle colour excitations is 

v»D^ J° + m 2 D v»F^ = Cun [J°] + Cfo v) , (25) 

where the linearised collision integral Cu n can now be evaluated explicitly, using 
the expressions for a and bj. The Yang-Mills equation remains unchanged at 
this order. For the collision integral, one finally obtains to linear order in the 
mean current, and to leading logarithmic order (LLO) in (lnl/g) -1 <C 1 

CW#](*,v) = 9f abc (al(x)5J^(x,v))\ A ^ LLQ 

= -1 I ^I(vy)J a °(xy) (26) 
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with the kernel 7(v,v') = <5 2 (v - v') - ^(v.v') 2 /-^ 1 - (v.v') 2 . Notice that 
the collision integral is local in coordinate space, but non-local in the angle 
variables. The rate 7 = ^-hxl/g is (twice) the hard gluon damping rate. 
In obtaining (p6|), we introduced an infra-red cut-off ~ gniD for the elsewise 
unscreened magnetic sector. This result has-been first obtained by Bodeker 
in ref. For alternative derivations, see ref. In addition, the source for 

stochastic noise £ in ( p5| ) can be identified as 

C„(z,v) = 9fabcal(x)SJ^(x,y)\ A=0 J=Q . (27) 

Making use of the basic correlator (||), its self-correlator follows to LLO as 

{C{x 1 v)C b {y,V))=2 1 mlTI{vy)5 ab 5{x-y) . (28) 

Notice that the strength of the correlator is determined by the kernel of the 
collision integral, e.g. by the dissipative process. Next, we show the close 
relationship to the fluctuation-dissipation theorem (FDT). 



9 Fluctuation-dissipation theorem 

It is well-known that dissipation in a quasi-stationary plasma (= no entropy 
production) is intimately linked to the fluctuations, a link which is given by the 
FDT. We shall argue that the FDT holds true for the above set of equations. 
To that end, following ref. El, we consider the coarse-grained kinetic equation 
and ask how the spectral functions of the noise source £ and of / have to 
be related to C[f] in order to satisfy the FDT. 

Within, classical transport theory, the FDT is implemented in a straightfor- 
ward wayt3. The pivotal element is the kinetic entropy S[f]. We consider small 
deviations of / from the equilibrium, / = f cq + Af. The entropy, stationary at 
equilibrium, is expanded to quadratic order in A/, S — S cq + AS. Defining the 
thermodynamical force in the usual way as F(z) = —5(AS)/5(Af(z)), where 
z = (x,p,Q), we obtain 

F(x,p,Q) = Af(x,p,Q)/f eq (29) 

for a classical plasma. The thermodynamical force fot-tiie quantum case follows 
from (||) using (|). Given C[f](z), it can be shownBlia that 

is the required noise-noise self-correlator for £ compatible with the FDT. For 
the particular example studied above, (|28| ) follows from inserting (p6| ) into 
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(|30|), proving that Bodeker's effective theory is compatible with the FDT. 
Furthermore, the correlator (A/A/)| t=t / can be derived along the same lines 
and agrees, to leading order, with (6 j '5 f)\t=t> derived from the Gibbs ensemble 
average. This guarantees that the formalism is consistent with FDT. 



10 Summary 

We have reviewed a semi-classical approach to derive effective transport equa- 
tions for QCD based on 'integrating-out' fluctuations about some mean fields. 
Most interestingly, the approach is applicable for in- and out-of-equilibrium 
situations, opening a door for future applications to out-of-equilibrium plas- 
mas. The formalism is consistent with the underlying non-Abelian mean field 
symmetry, and allows for systematic expansions controlled by a small plasma 
parameter which is at the basis for any reliable computation. 

For thermal plasmas we have discussed how the collision integral as well 
as the necessary noise source follow explicitly from the microscopic theory 
to leading logarithmic accuracy, reproducing Bodeker's effective theory, and 
establishing the consistency with the fluctuation-dissipation theorem. 

As a final comment we point out that the effective theory is the same for 
a classical or a quantum plasma, differing only in the equilibrium distribution 
function, and hence in the corresponding value for the Debye mass. The sole 
'quantum' effect which entered the computation resides in the non-classical 
statistics of the particles, which is all that is needed to correctly describe a 
hot quantum non-Abelian plasma close to equilibrium at the present order of 
accuracy. 
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